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We discuss how multichannel Kondo physics can arise in the setting of a localized level coupled 
to several bosonic Tomonaga-Luttinger liquid leads. We propose one physical realization involving 
ultracold bosonic atoms coupled to an atomic quantum dot, and a second, based on superconduct- 
ing nanowires coupled to a Cooper-pair box. The corresponding zero- temperature phase diagram 
is determined via an interplay between Kondo-type phenomena arising from the dot and the conse- 
quences of direct inter- lead hopping, which can suppress the Kondo effect. We demonstrate that the 
multichannel Kondo state is stable over a wide range of parameters. We establish the existence of 
two nontrivial phase transitions, involving a competition between Kondo screening at the dot and 
strong correlations either within or between the leads (which respectively promote local number- 
and phase-pinning). These transitions coalesce at a self-dual multicritical point. 
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Introduction. Magnetic impurities can drastically al- 
ter the low-temperature properties of metals, leading to 
anomalous temperature dependence in, e.g., the heat ca- 
pacity, resistance, and magnetoresistance. These prop- 
erties, collectively known as the Kondo effect [1 , are 
exhibited when the impurity interacts antiferromagneti- 
cally with the conduction electrons of the metal, and are 
due to the dynamic screening, as T ^ 0, of the spins of 
the individual impurities by a cloud of electrons. The 
low-energy scattering properties of each screened impu- 
rity are those of a bound but spin-polarizable impurity- 
cloud singlet^ whose effects on the electron gas can be 
described using Fermi-liquid theory [2 . For a ferromag- 
netic impurity-cloud interaction, on the other hand, the 
spins of the impurity and its polarization cloud align in 
one of two degenerate triplet ground states. A quantum 
phase transition separates the ferro- and antiferromag- 
netic (i.e., Kondo) cases. 

A striking generalization of this "1-channel" Kondo ef- 
fect is the multichannel version [3 , in which each im- 
purity couples separately to conduction electrons that 
propagate in N{> 1) channels (e.g., distinguished by or- 
bital angular momentum). When N > 2s (with s being 
the spin of the impurity), the conduction electrons over- 
screen the impurity, and the low-energy behavior can no 
longer be described by Fermi-liquid theory; instead, the 
thermodynamic properties follow anomalous power-laws 
governed by a quantum critical point [4^ . In metallic sys- 
tems, observing the A/'-channel Kondo effect has proven 
demanding [5]. Although advances in nanoscience have 
enabled the exploration of Kondo physics in the more 
controlled setting of 2DEG-based quantum dots con- 
nected to leads (which realize the channels), even here 
the engineering of multichannel Kondo phenomena re- 
mains a challenge. This is primarily due to the dif- 
ficulty in preventing interchannel hybridization, which 
gives rise at low energies to a single composite channel 



that screens the impurity via the 1-channel Kondo effect. 
For instance, signatures of the two-channel Kondo effect 
have recently been observed in a quantum-dot-based set- 
ting [7^, but required fine tuning to prevent hybridization. 

In view of these challenges, it is desirable to explore 
the TV-channel Kondo effect and its onset in a more read- 
ily controllable setting: this can be achieved, e.g., using 
leads having tunable interparticle correlations. For the 
1-channel Kondo effect the possibility of such control was 
demonstrated in Ref. [8 in the context of a quantum dot 
coupled to a Tomonaga-Luttinger liquid (TLL) lead [9]. 
The strength of the repulsive interactions in the lead (as 
encoded in the TLL parameter K) was found to tune 
the position of the Kondo-to-ferromagnetic phase tran- 
sition. Motivated in part by this result, in this paper 
we explore the case of N TLL leads coupled to a quan- 
tum dot; this case is expected to exhibit richer physics 
arising from the possibility of not only intra-lead but 
also inter-lead correlations. We focus on the case of 
bosonic leads because for them a wide range of K val- 
ues can be experimentally accessed with relative ease: 
e.g., ultracold bosons with short-range repulsions have 
K > 1 [10 , whereas ones with dipolar interactions have 
K < 1 pj. We suggest two concrete realizations: one 
uses cold atoms [121 US] ; the other uses superconducting 
nanowires and a Cooper-pair box [14 . Systems of ul- 
tracold atoms are especially well suited to the study of 
multichannel Kondo physics because their interactions 
are highly controllable and tunable [10 , and extraneous 
noise can be mitigated. 

Our main theoretical results concern the competition 
between the A^-channel Kondo effect and interactions in 
the leads; these interactions either suppress lead-dot tun- 
neling or generate inter-lead phase-locking (which would 
short-circuit the dot). Our analysis yield a phase diagram 
containing four distinct phases (see Fig.jsj below) and ex- 
hibiting a pair of unusual phase boundaries, which meet 
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at a self-dual multicritical point. As discussed below, an 
important advantage of our proposed experimental real- 
ization of the A^-channel Kondo effect is that it is robust 
against interlead hybridization. We discuss the condi- 
tions for observing the A/'-channel Kondo effect using cold 
atoms, and suggest that the unusual phase boundaries 
may be accessed experimentally using dipolar bosons. 
Kondo and resonant-level models. We consider the gen- 
eral anisotropic Kondo Hamiltonian Hk = ^leads + ^int 7 
where i^ieads = Xl/c k c^^ describes a free-electron 
gas (or a Fermi liquid) and i^int describes the coupling 
to the impurity spin, which is located at r = 0: 



^int = 



r 



^+4(o)c^(o) + h.c. ^4(oK(o), 



where Ccr(O) annihilates a conduction electron of spin cr at 
the impurity location, the Pauli matrices S act on the im- 
purity spin state, and and J^, respectively, the ampli- 
tudes for the lead-dot Ising and spin-flip processes. The 
antiferromagnetic case, > 0, leads to Kondo screen- 
ing; for < and \J±_\ < \Jz\, the impurity couples 
ferromagnetically to the conduction electrons. 

We shall primarily be concerned with a model that 
is equivalent to the Kondo model, viz., the interacting 
resonance-level model (iRLM). This model consists of a 
localized level d (i.e., a quantum dot) at the Fermi energy 
hybridized with N channels of spinless noninteracting 
conduction electrons c^j^, together with a short-ranged 
repulsion between dot and lead electrons: i^iRLM = 
Pleads + ^onsite + ^int- Here, the Icads are described 
by i^ieads = Z^kZln^i ^(k) 4k ^nk' ou-dot potential 
by i^onsite = B d^ d^ and the dot-lead couplings by 

Hr^r.t=V{d^c„iO)+h.c.) +U{dU- 1 ) (4(0) c„(0) - 1 ) (1) 

where V ^ J± and U ^ Jz -\- const. The canonical 
transformation that maps the Kondo model on to the 
iRLM consists of identifying the spin density waves of 
the Kondo model with the particle density waves of the 
iRLM. The analogy is most evident as /7 ^ oo: the pres- 
ence of an electron on the dot ensures the absence of 
electrons near the dot, and vice versa, i.e., an anticor- 
related state resembling the Kondo singlet. The Kondo 
effect manifests itself in the iRLM via an enhancement, 
as T ^ 0, of number fluctuations on the dot. 

We now turn to the iRLM with TLL leads, governed by 
the Hamiltonian i^ieads = Z]fe^(^)4^/c' which {bk} 
are bosonic fields that describe free collective phonon 
modes that disperse linearly [9]: E{k) ~ \k\. As discussed 
in Ref. [8j, this version of the iRLM can also be mapped 
on to the standard Kondo model, but the location of 
the ferro- to antiferromagnetic transition depends on K. 
The boundary between the ferro- and anti-ferromagnetic 
phases occurs at = hVs{^/l/2K — 1), in which Vs is the 
speed of sound in the leads. 

In the A/'-channel generalization of the iRLM, each 
channel couples independently to the impurity via H^^^. 



By inverting the iRLM-Kondo mapping, one arrives at 
the A/'-channel Kondo Hamiltonian studied in Refs. OH], 
in which each of A' channels couples independently to 
the spin. However, the interchannel particle transfer 
terms (i.e., hybridization) ilp"' ~ ^tai^) ^c^P ^n' [3^^)^ 
which destabilize the A/'-channel Kondo effect, do not 
arise in the iRLM. This is because the particle-number 
(or "charge" ) sector of the leads in the iRLM — the model 
of physical relevance here — maps on to the spin sector of 
the equivalent Kondo model. In contrast, the particle 
sector of the leads in the equivalent Kondo model has 
no physical significance in the iRLM and, accordingly, 
cannot couple to the dot. 

Realizations of the iRLM. Our first realization involves 
ultracold atoms, and extends the ideas of Ref. [13 . A 
star-shaped pattern of A' one-dimensional leads meet- 
ing at a point can be constructed by passing a laser 
beam through a phase mask or spatial light modulator 
(SLM) [15 . Such a device is a sheet of glass of spatially 
varying thickness, which distorts flat wavefronts, giving 
rise to a prescribed intensity pattern at a "screen" some 
fixed distance away. Algorithms for the construction of 
appropriate phase patterns are discussed in Refs. [151 [16]. 




FIG. 1: Candidate cold-atom setup. The spatial light 
modulator (SLM) distorts the wavefronts of laser A so as to 
create a star-shaped pattern at the "screen'.' Atoms are con- 
fined to the screen using two lasers (B) propagating at a rel- 
ative angle ^ ^ tt so as to create an optical lattice of spacing 
L = A/ [2 sin(^/2)]. Lasers C and D are used to confine atoms 
in hyperfine state b at the "dot" (dark blue region). 

As for the dot itself, it can be realized as follows [13]. 
Suppose that the atoms discussed in the preceding para- 
graph are in a hyperfine state a. A tightly-confining trap 
for a different hyperfine state b is now created at the in- 
tersection of the leads. An atom in state a can make a 
Raman transition to state 6, and vice versa; while in 6, 
it is confined at the "dot'.' The Raman transition thus 
creates a lead-dot hopping amplitude. This setup real- 
izes an iRLM having the following couplings: K is de- 
termined by the scattering length gaa for atomic state a, 
U by the a ^ b scattering length Qah^ and V by the am- 
plitude (i.e., effective Rabi frequency) A of the Raman 
transition. Double occupancy of the dot is prevented by 
a large repulsive interaction ^55 between atoms on the 
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dot (with gaaidah)- All interactions are tunable 

via a Feshbach resonance. The direct interlead hopping 
amplitude is governed by the intensity of the laser that 
traps a-state atoms at the intersection of the leads. 

One can realize a similar model in a mesoscopic set- 
ting by using superconducting nanowires as the leads, 
together with a Cooper-pair box [14 — i.e., a supercon- 
ducting island that holds at most one Cooper pair — as 
the dot. The normal modes of the leads are plasmon 
excitations, and the Hamiltonian for the box- lead sys- 
tem has the same form as that for the cold atom system, 
provided the leads are connected to the dot via Joseph- 
son couplings (which determine V). The coupling U is 
determined by the lead-dot Coulomb repulsion. 

Analysis of the model. In addition to processes that 
involve the dot, we account for those in which bosons 
hop directly between leads. It is useful to write 
the Hamiltonian for the uncoupled leads as i^ieads = 

ElLi !o dx[K{dJr.? + K-\d,An)% where the den- 
sity fluctuation modes of the TLLs are given by the 
operator pn{x) ~ dx(\>n{p^^l^ and the canonically con- 
jugate momenta by dxOi{x). Direct hopping processes 
between the leads can be described using the boson 
annhilation/creation operators tpn{x) ~ e'^^''^^^\x=o at 
the end-points of the semi-infinite TLLs: i^tunn = 
t En,n'(e^^^"^°^"^"'^°^^ + h-c-)- At low frequencies, the 
ground state of the complete system consists of either 
uncorrelated wires (i.e., the disconnected fixed point, 
DFP, t = 0), OT N maximally correlated wires (i.e., 
the connected fixed point, CFP, t ^ oo). The CFP 
manifests itself via the mutual pinning of the phase 
fields at the junction: On{x^t) = On'{x^t)\x=o for all 
pairs {n^n') [17 . Additionally, current conservation de- 
mands that St^n(0, t) = 0. In the language of 
the renormalization group (RG), the hopping t fiows 
to (i.e., DFP) for < 1 ^8j, and to strong cou- 
pling (CFP) for > 1 [m [H]. An RG analysis 
around the CFP for the backscattering amplitude A gives 
d\/dl = {1 - (2{N - l)K/N)} A. Therefore, the CFP is 
stable for K > N/ (2(7V - 1)). Both the DFP and the 
CFP are stable for N/[2{N - 1)] < < 1; for K in this 
interval there must therefore be a quantum phase tran- 
sition at some / 0, separating the ground states of 
N uncorrelated and N maximally correlated wires. This 
transition is analogous to that exhibited by a quantum 
Brownian particle on an [N — l)-dimensional triangular 
lattice in the presence of Ohmic dissipation [19]. 

The novel feature of the A/'-wire phase diagram shown 
in Fig. [2^ is that the CFP is stable against weak backscat- 
tering close to the junction, even for < 1 (in contrast 
with the two- wire case ^). This occurs because the 
boundary conditions on the fields at the junction imply 
that a wavepacket arriving at the junction from any one 
wire meets an effective composite TLL, comprising the 

- 1 other leads, for which i^eff = {N - l)K > 1. 



a. b. 




FIG. 2: (a) Phase diagram for N leads with interlead hopping 
t and TLL parameter K. The phase boundary (thick line) 
separates the disconnected (i.e., DFP, t = 0) and maximally 
connected (i.e., CFP, t oo) fixed points, and has a self-dual 
point (M). (b) Kosterlitz-Thouless flow for lead-dot couplings 
at fixed (K, t) . The left separatrix (thick line) demarcates the 
boundary between the ferromagnetic and Kondo phases. 



The locking of the phase field in the first wire to that 
of the composite TLL precludes any chemical-potential 
drop across the junction, although only a fraction of the 
incoming current enters any individual lead of the com- 
posite TLL [17 . Said another way, the local inertia of 
the phase fields strongly suppresses backscattering events 
involving high- momentum phase fiuctuations. This phe- 
nomenon is dual to the enhanced inertia in the number 
fields at the endpoints of the wires at the DFP ioi K < 1, 
which results in a power-law suppression of the tunneling 
density of states (TDOS). 

We now incorporate the dot by coupling it to the N 
leads via H^!^^ [see Eq. ([T])], and develop the phase dia- 
gram via an RG analysis around the CFP and DFP. Near 
the CFP, we find the following scaling equations: 

to second order in all couplings, where the couplings 
have been scaled by the high-frequency cutoff uJc = Vs/£, 
(where ^ is the healing length |13j). In addition, we 
find that dX/dl = {1 - {2{N - 1)K/N)} X-}- {J'^ /cjc). By 
shifting to = + 1 - 2{N - l)K/N, Eqs. ^ as- 
sume the well-known Kosterlitz-Thouless form [20j (see 
Fig. ^) , with a Kondo temperature scale given by Tk ~ 
cjcG Z'^^. Thus, Jz is found to be RG-marginal at first 
order but RG-relevant at second order, and independent 
of K and N. On the other hand, J± has the same scal- 
ing dimension as A, and is thus dependent on K and N. 
For K > [N/2{N - 1)], even though J±_ is RG-irrelevant 
(from its scaling dimension), it can turn relevant, due 
to the growth of J^. The scaling equation for J± ad- 
mits a nontrivial fixed point at Jz = [2{N — 1)K/N] — 1. 
For Jz > Jz, all fiows lead to the A^-channel Kondo fixed 
point; for Jz < Jz, fiows lead to zero dot-to-lead hopping. 

IfK < [N/2{N - 1)], the RG analysis about the CFP 
is invalid; one must instead analyze the Kondo couplings 
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around the DFP. The scaling relations near the DFP are 

dJj ^ 2 = {i_ 1_ 

dl ~ dl ~ \ K 



(3) 



As with A, the flow for t acquires a positive contribution 
of order J\/uOc from the dot-mediated hopping (i.e., the 
dot promotes interlead hopping, as one might expect). 
By shifting to = + 1 - (1/i^), Eq. ^ assumes 
the Kosterhtz-Thouless form. Similarly, for < 1 the 
scaling equation for has a nontrivial fixed point at 
Jl = (1 - K)/K. For > j;, ah flows lead to the 
A/'-channel Kondo fixed point; for Jz < J^^ flows lead to 
zero dot-to-lead hopping. 
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FIG. 3: Phase diagram of the dot-lead system. The discontin- 
uous critical surface separates the Kondo and ferromagnetic 
phases. The curved vertical ribbon (orange) is the phase 
boundary between the DFP and CFP phases (see Fig. |2|. 
The discontinuity of the critical surface shrinks to zero at the 
multicritical point M. The diagram has four phases — I: com- 
pletely decoupled wires; II: wires coupled only through dot; 
III: wires connected both directly and through dot; IV: wires 
connected but decoupled from dot. The transitions between 
them (shown in insets A, B, and C) are described in the text. 

Bringing together the flows of (t^J^^Jz) yields the 
three-dimensional phase diagram shown in Fig. [3j The 
tuning of t and/or K allows one to access two non- 
trivial transitions between phases that have opposing 
characters in hoth their Kondo coupling to the dot and 
their direct interlead hopping. One is a transition be- 
tween phase II (in which A^-channel Kondo physics dom- 
inates TDOS suppression) and phase IV (in which Kondo 
physics is suppressed by local phase pinning); see Fig.|3^. 
The other is a transition between phase I (in which 
Kondo physics is dominated by TDOS suppression) and 
phase III (in which Kondo screening overcomes local 
phase pinning); see Fig. [3p. In addition, we find a mul- 
ticritical point (see point M in Fig. |3|^) at intermedi- 
ate coupling in t; this occurs when = J* [so that 




FIG. 4: RG flows in the t- J± plane for Nl2{N - 1) < K < 
^N/2{N-1). Here, < |J. - J:|,|J. - J.|. I. J. < 
Jz,Jz- hopping via the dot is irrelevant, and the dot stays 
disconnected. II. Jz < Jz < Jz' for J± 0? hopping via the 
dot is irrelevant for t < t* and relevant otherwise. However, 
sufficiently large J± can drive t past t* toward the regime 
in which both t and J± grow at low energies. HI. Jz > Jz- 
hopping via the dot is relevant on both sides. Figures I and 
HI would remain identical for ^N/2{N - 1) <K <l, while 
II would become inverted about the dashed line. 



K = y]v72(]V^^]. Point M coincides with the self- 
dual point [19 in the phase boundary of intermediate- 
coupling fixed points (see Fig. |2^) , and involves a com- 
promise between the competing tendencies of TDOS sup- 
pression and local phase pinning. For the special case of 
N = 2^ the phase boundary in Fig. [3^ becomes a marginal 
line at = 1 (i.e., the Tonks-Girardeau gas [H]), and 
the point M becomes a multicritical line at = J* = 0. 

The four phases are characterized by the following bare 
two-lead transmission coefficients across the junction: 
G2 = (i.e., minimal) at the DFP and G2 = 4:K/N^ (i.e., 
maximal) at the CFP. These bare coefficients acquire 
power-law corrections of order (DFP) and 

(CFP), arising from direct interlead scattering [18 ; here, 
T represents an energy scale (e.g., temperature) that cuts 
off the RG flows, and {u^ ja) are exponents determined 
by the leading irrelevant perturbations around the cor- 
responding fixed point. In experiments with cold atoms, 
these power-law contributions should be detectable via 
real-time lead dynamics [17 . 

Number fluctuations on the dot can be accessed via a 
non-destructive measurement scheme such as, e.g., that 
suggested in Ref. [21] for the Bose-Hubbard model: in 
such a scheme the dot would be located in the waist 
of a high-finesse optical cavity that has a resonance fre- 
quency near an optical transition of the hyperfine state 
b (but not of a). A fixed-number state on the dot (i.e., 
the unscreened spin) would merely shift the cavity's res- 
onance; by contrast, a fluctuating- number state (i.e., the 
Kondo state) would lead to a double-peak structure in 
the transmission spectrum of the cavity, with peaks cor- 
responding to an empty dot and to an occupied dot. Ac- 
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